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The only information available about an alleged source of entangled quantum states is the amount 
S by which the Clauser-Horne-Shimony-Holt (CHSH) inequality is violated: nothing is known about 
the nature of the system or the measurements that are performed. We discuss how the quality of 
the source can be assessed in this black-box scenario, as compared to an ideal source that would 
produce maximally entangled states (more precisely, any state for which S = 2\/2). To this end, we 
introduce several inequivalent notions of fidelity, each one related to the use one can make of the 
source after having assessed it; and we derive quantitative bounds for each of them in terms of the 
violation S. We also derive a lower bound on the entanglement of the source as a function of S only. 

PACS numbers: 03.67.-a, 03.65.Wj, 03.65.Ud 



I. INTRODUCTION 

A device, allegedly generating pairs of entangled par- 
ticles, is for sale. Obviously, the potential user wants to 
check that entanglement is indeed being generated be- 
fore buying it; but just as obviously, the vendor does not 
want to open the device and reveal its fabrication. For 
classical devices, such a situation would lead to a com- 
plete impasse. Not so, however, for quantum devices: 
Bell's inequalities can act as entanglement witnesses ir- 
respective of the nature of the system under study or 
of the kind of measurements that are being performed. 
Thus suppose that the vendor provides the user with two 
additional boxes, the measurement devices. Once more 
the vendor does not want to open the device and reveal 
its fabrication. Suppose in addition that the user can 
choose the measurements: the measurement devices have 
a knob whose positions correspond to allegedly different 
measurements (Fig. [TJ. By operating these devices, the 
user can reconstruct the statistics P(a, b\A, B) of the ob- 
served outputs a and b, conditioned on each choice of 
knob positions A and B. If the statistics violate some 
Bell inequality, and the measurement has been performed 
in such a way as to avoid signaling between the measure- 
ment boxes, then the user is convinced that the source is 
indeed producing entangled pair. 

The possibility of such an assessment is already re- 
markable. However, the user cannot be satisfied with 
knowing that there is "some entanglement" : what is 
needed is a quantitative estimate on how good the source 
actually is. The amount of violation of a Bell's inequality 
can provide such a quantitative criterion, provided it is 
translated into the meaningful figure of merit: fidelity or 
trace distance to the ideal state, or some entanglement 
measure... The goal of this paper is to provide such quan- 
titative estimates, when the Bell inequality under study 
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FIG. 1: Device-independent state estimation: the quality of 
an unknown source of entangled pairs should be established 
using unknown measurement devices. The only available in- 
formation are the statistics P(a, b\A, B) of the outcomes (a, 6) 
for measurement settings A, B. The figure represents the par- 
ticular case studied in this paper, where both the choice of 
measurement settings and the outputs are binary. 



is the CHSH inequality derived by Clauser, Home, Shi- 
mony and Holt P|. 

This work is inspired by "device-independent quantum 
key distribution" 0,0], in which the amount of violation 
of the CHSH inequality is used to bound the information 
of an eavesdropper without making any hypothesis on 
the internal workings of the devices. It is also related to 
the concept of "dimension witness": sufficient violation 
of some Bell inequalities can guarantee that the quantum 
state has a minimum dimension @, 0, Q . One possible 
application of the present work could be to devise im- 
proved self testing of quantum computers 0, [1] • 



II. FORMULATION OF THE PROBLEM 
A. Ideal states 

As we said, we restrict to the case where the user ap- 
plies only two measurement settings on each particle and 
the outcome is binary. In this case, there is only one 
Bell inequality, namely CHSH [§]. We further restrict 
our study in considering only the observed violation S b s 
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of CHSH as quantitative measure, being aware that the 
statistics P(a, b\A, B) contain further information that 
might improve the estimates. 

Since the source will be characterized by a single scalar 
quantity, the set of ideal states is the set of states 3? 
such that S = 2^/2 is achievable. This set has been 
fully characterized [l(| EH : it consists of all pure states 
of the form X)j c il^'j) where is a two-qubit max- 
imally entangled state in a four-dimensional subspacc, 
i.e., |%) = ^(|2j - 1, 2j - 1) + \2j, 2. 7 )) up to local uni- 

taries. Since the relative phases of the Cj do not play 
any role in the violation, we must add mixed states to 
the set. It is easy to verify that the most general such 
state can be written as $ = UaUb^ + ® a ^ a^b where 
|$ + ) = -i=(|00) + 1 1 1 ) ) is a two-qubit maximally entan- 
gled state, <r is an arbitrary state, and Ua, Ub are ar- 
bitrary local unitarics. In their work on device-testing, 
Mayers and Yao (MY) [t| chose their reference states as 
those that could be written in the above form with a 
pure, i.e., even though they did not refer to Bell inequal- 
ities, they where considering all pure states that violate 
CHSH maximally. 



and one is then led to search for bounds of the form 

F M Y{p)>J r MY{S ohs ). (5) 

For Tmy we obtain tight lower bounds if the state is re- 
stricted to consist of two qubits, or (modulo a conjecture 
of Gisin and Peres) if the state is restricted to be pure. 
Putting such hypotheses on the source goes against the 
philosophy of the black box scenario, but it allows us to 
get a mathematical grasp of the problem. When no re- 
strictions are put on the state, we do not even have a 
lower bound on Tmy ■ However it is possible to intro- 
duce other notions of fidelity (see below) which have a 
clear operational meaning, and for which lower bounds 
can be computed without any hypothesis on the source. 

Yet an alternative approach to the source characteriza- 
tion problem would consist in looking for a lower bound 
to the entanglement of the state p: 

E{p) > S(S ohs ) , (6) 

where E is an entanglement measure, such as the entan- 
glement of formation, of distillation, etc... [HI]. Below 
we obtain lower bounds on £. 



B. Figures of merit 

The distance between the actual source state, with 
density matrix p, and the closest ideal state 3>, is conve- 
niently measured by the trace distance [H, EH : 

Smy(p) = min£(p,#) , (1) 

where 8(p, = |Tr|p — The trace distance has 
a clear operational interpretation: in whatever task, p 
will behave differently from $ with probability at most 
S(p, <I>). In other words, the real source will differ from 
an ideal source with probability at most S(p). 

The problem we have set out to solve is thus to find a 
bound of the form 

6my(p) <V M Y(S ohs ). (2) 

This bound can in principle be obtained by solving the 
following optimization problem: 

V my (Sobs) = max |min5(p, *)i , (3) 

p:S ma , x (p)>S obs I- * > 

where 5 , max (p) is the maximum CHSH violation that can 
be obtained by measuring state p. 

Deriving lower bounds, let alone tight lower bounds, 
for T>my turns out to be much harder than we initially 
anticipated. In practice, it is simpler to work with the 
Fidelity rather than the trace distance, the two mea- 
sures being related by S < \/l — F [l3| ■ In analogy with 
cq. |T]), we define 

Fmy (p) = max F(p, *) = max f Try/ p 1 / 2 ^>p 1 / 2 ) ^ (4) 



C. Warm-up: solution assuming two qubits 

As a nontrivial warm up exercise, let us compute the 
bound cq. ([5|) under the assumption that the source 
emits a pair of qubits and that the measurements are 
von Neumann measurements. This is an undue restric- 
tion for the black-box scenario; we present this calcula- 
tion because its result is interesting in itself, and will be 
an important tool for the main discussion. 

In this case the set of ideal states is well known: only 
the maximally entangled states 3? = UaUb^ + U a II^ 
violate CHSH maximally. Therefore, Fmy(p) = 
max$ F(p, «fr) reduces to the so called singlet fidelity of 
p. Our approach consists in fixing the singlet fidelity 
of p, and computing S max (p). To this end we use the 
spectral decomposition of the Bell operator 

B = (A + A') <g> B + (A - A') g> B' . (7) 

First note that if F(p) < ^ the state cannot be entan- 
gled, CHSH cannot be violated, and the bound S max = 2 
can be trivially achieved by the degenerate measurement 
A = A' = B — B' — 1. If the inequality is violated, the 
operators A, A , B, and B' must be linear combinations 
of the three Pauli matrices. Then the spectral decomposi- 
tion B = J2i )>i\$i)(®i\ has the following properties [l5|: 
the |<!>i) are a Bell basis (i.e., a basis of maximally entan- 
gled states) and the eigenvalues are {Ai, A2, — A2, — Ai} 
with Tr(£? 2 ) = 16, i.e., 

A? + \ 2 2 = 8 , (8) 

which implies the Cirelson bound |Ai| < 2^/2 [lj|. 
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Therefore, for a given B we have 

S(p) = Tr(pB)=J2M$i\p\*i)- (9) 

i 

Suppose for dcfinitcncss Ai > A2 > 0. Then, keep- 
ing F(p) fixed, S(p) is maximized by choosing |$i) 
such that F(p,&±) = F(p). Whereupon we have that 
S m ax(p) < XiF(p) + A 2 (l — F(p)) because the two other 
eigenvalues are non-positive. Using eq. I|8|). we can set 
Ai = 2\/2cosa; and A 2 = 2y/2 sinx. The well-known 
bound max x (fl cosi + b sin x) = \/a 2 + b 2 then leads to 
SmM < 2 % /2^F(p) 2 + [l-F(p)] 2 . Finally (Fig.©: 

Fmy(p) > (l + VIW2] 2 - l) /2 (qubits). (10) 

This bound is tight, being achieved by pure non- 
maximally entangled states \ip) = cos#|00) + sin#|ll). 

Indeed, for these states 5 max = 2^1 + sin 2 (26>) [H El 

and the singlet fidelity is F = \ (V>|$ + ) | 2 = |(1 + sin(20)). 
Furthermore for pairs of pure states, we have the strict 
equality 5 = vl — F, hence eq. (fTQ|) leads to a tight 
bound for the trace distance as well. 



more careful, because the above construction does not 
imply: 

Fmy(p) > ^p a pF{p afj ) (probably wrong), (13) 

where F{p a p) is the singlet fidelity of p a p. The reason 
is that in the MY approach, the state must be brought 
close to a reference state using local unitary operations 
V A <£> Ub- Let U a be the restriction of U a to the 2x2 
block indexed by a; and similarly for Up; and let 3> a /3 be 
the maximally entangled state of two qubit such that 
F(& a/3 ,p al3 ) = F{p af} ) is the singlet fidelity of p af3 . 
Now, there is no guarantee that Ua and Ub exist, such 
that U a <X> Up&cpU^ ® = 3> + for all a and f3, as is 
required to obtain a reference state according to the MY 
definition. Moreover, the MY definition of fidelity is a 
comparison with the whole state 3> + ® rr, not only with 
the two-qubit component <I> + . In order to make sense of 
eq. (fl3)) we will introduce different definitions of fidelity 
below. Before turning to that, we present the case of 
pure states of arbitrary dimensions, for which the MY 
fidelity can be computed. 



III. BOUNDS ON THE FIDELITY TO THE 
CLOSEST REFERENCE STATE 

A. Structure of the Bell operator 

For any two dichotomic operators A and A', one can 
find a basis such that both operators are block-diagonal, 
where each block is a 2 x 2 matrix (see e.g. [||). So one 
has A = Y, a A a and A' = J2 a A' a where A a = Tl a An a , 
A' a = H a A'H a and II Q are orthogonal projectors onto 2 
dimensional spaces. Of course, a similar decomposition 
holds for Bob's operators. Therefore, the Bell-CHSH op- 
erator can be written as 



<*,0 ■ 



(11) 



a, 



where B a _p = J2i X^\^^)(^^\ are orthogonal two- 
qubit operators with the same properties as above. 
Therefore 

S(p) = ^2p a pTr(p a/3 B af )) = ^2p a pS(p af) ) , (12) 

a,0 a,/3 

where PapPap = ® n^pn a ® lip and p a a is a normal- 
ized two-qubit state. 



B. A complex problem 

Given (p~2|) , it may seem that the extension of our result 
to arbitrary dimensions is just a matter of convex opti- 
mization. A closer look shows that one must be much 



C. Solution under the restriction to pure states 

Let us assume that we know that the source emits 
a pure state (again an undue restriction for the black 
box scenario). Using the Schmidt decomposition I*) = 
^ fc Afc|fc,fc) with the Schmidt coefficients in decreasing 
order > A^+i > 0, any pure state can be rewrit- 
ten as = ^ ^/pj(cj\2j, 2j) + Sj \2j + l,2j + 1» with 
c 2 + s 2 = 1. The MY fidelity can be computed exactly 
(Appendix [A"|) . and one finds 



F, 



AIY 



(*) = £ 



(\j + Xj+iY 



&J <a±2l!. (i4) 



This should now be related to S'max(^')- For states of ar- 
bitrary dimension, there is no known analytical expres- 
sion for the maximal violation of CHSH. However, for 
pure states ^ there is a long-standing conjecture by Gisin 
and Peres [l9j . whose validity has never been disproved 
by numerical checks [2(j. According to this conjecture, 
the ordered Schmidt decomposition defines the natural 
block-structure of the CHSH operator. This implies 



C 3 S J 



(15) 



Combining this conjecture with eq. (|14[) we find that for 
pure states the accessible points in the (Fmy , •S'max) plane 
are convex combinations of points on the curve given by 
equality in eq. (TIT))) , yielding (Fig. [2]): 



Fmy(*) > [S ohs + 2V2 - 4] 

(pure states, modulo Gisin-Pcrcs conjecture). 



(16) 
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This bound is tight if wc allow the dimension d to become 
arbitrarily large (otherwise, the ordering of the implies 
constraints on the possible values of {pj,Cj, Sj}). More- 
over, this bound is weaker than the one obtained under 
the assumption of two-qubits. Though not astonishing 
in itself, this feature is new: in device-independent quan- 
tum key distribution, the bound for collective attacks is 
already optimal in the two-qubit case @, [|[ . 



-S 0.8 



o 0.6 



; 0.4 



S 0.2 



MY qubits 






V- 








/"MY pure = 


LOCC 








! 2.2 


2.4 2.6 


2.8 



FIG. 2: (Color online) Lower bounds on the fidelity as a func- 
tion of the observed violation S'obs of the CHSH inequality. 
From top to bottom: Fmy assuming two qubits (|10|l ; Fmy as- 
suming pure states and the Gisin-Peres coniecture (|16[l . equal 
to Flqcc (HHI; and F LO 0). 



D. Black-box bounds for other fidelities 

The MY fidelity is defined to suit the black-box sce- 
nario. However, other definitions of fidelity may be mean- 
ingful. Here, we consider fidelities defined as 



F c (p) = maxF(A(p),*+). 



(17) 



where £ is a set of completely positive maps which map 
the Hilbert space of p onto a 2 x 2 dimensional Hilbcrt 
space, and which cannot increase the entanglement. Fc 
can be thought of as the best singlet fidelity obtainable 
under single shot purification of p to a two-qubit en- 
tangled state, using only operations that belong to the 
family C. Wc will consider the case where C consists of 
all the completely positive maps that can be realised by 
Local Operations (C — LO) or by Local Operations and 
Classical Communication (C = LOCC). 

These new notions of fidelity shed a different light on 
the task of source characterization. Indeed the Mayers- 
Yao fidelity and trace distance compare the state pro- 
duced by the source to the closest ideal state, thereby 
establishing how much the real state and the idea state 
would differ in applications. The new fidlities Fc are 
relevant to another scenario in which the user may try 
to improve the source by acting locally on the two sub- 
systems, for instance by opening the boxes containing 
the measurement devices and tinkering inside them. But 
before buying the source, the user wants to perform a 
fast black-box check to ascertain what will be the perfor- 
mance of the improved source. In other words, by mea- 
suring S'obs in a black-box scenario, the user can assess 



how well the source would perform in other scenarios. In 
this sense, the bounds we derive for these new fidelities 
are real black box statement, which do not make any hy- 
pothesis on the state p and on the measurement devices. 
These new fidelities are related by 

Fmy < Flo < Flocc ■ (18) 

In Appendix IB! we prove that Fmy(^) = Flo(*&) for 
pure states while for mixed states there are explicit cases 
of strict inequality. We also show, see Fig. ^ that: 

Flo(p) > - i - [S* obs - 2] , (19) 



2(\/2-l) 



LOCC 



(P) 



> 



1 



Sobs + 2V2-4 



(20) 



4(V2-1) 

The bound ([!§)) on Flo is obtained by exhibiting an 
explicit LO strategy. The proof is lengthy and we give 
it in Appendix [Cj We just note here that this bound 
is surely not tight, since it reaches the over-pessimistic 
F = for S = 2. 

The bound (|2"D|) for Flocc is the same one obtained 
for Fmy on pure states, eq. (|16[) : we do not know whether 
this bound is tight. The proof goes as follows. The de- 
composition eq. (fTTj) of the Bell-CHSH operator gives 
us a natural method for projecting p onto a 2 qubit 
space: particle A is projected in the LT a spaces and par- 
ticle B in the IL3 spaces. Using CC, the actual block 
(a, (3) is made known in both locations. The result 
of this completely positive map is a state with fidelity 
F = Pa/3 J2 a f3 F{p a p), i.e., we obtain a convex combina- 
tion of points on the curve eq. (|10[) . The bound (|2TH) is 
recovered by noticing that, in the LOCC scenario, all the 
blocks (a, (3) for which S(p aj3 ) — 2 can be brought to 
have F = 1/2: indeed, for the blocks where they observe 
5 = 2, Alice and Bob can swap their local states with 
those of ancillas prepared in a pure product state. 

IV. OTHER FIGURES OF MERIT 

The core of our work involved using the fidelity as a 
figure of merit. Here, we present the consequences of the 
bounds obtained on the fidelity for other figures of merit. 

A. Relation with trace distance 

Even if fidelity bounds were found to be tight, the 
tightness of the bound 8 < \/\ — F on the trace distance 
would follow only if the states that saturate the bound are 
pure. However, we are already able to conclude that the 
bounds I?(S' bs) for the trace distance 8 put very strin- 
gent constraints on the quality of the source. 

For instance, our strongest bound eq. (fTU|) leads to a 
tight 8 = VI - F. If we insert S'obs = 0.99 x 2\/2, we 
obtain <5 sw 10%. If the user requests the error rate to 
be below 1%, the vendor will have to produce extremely 
good sources — better than any currently available one. 
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B. Relation with entanglement measures 

The bounds for all the T c also provide lower bounds 
on the entanglement of p. Indeed, consider any entan- 
glement measure E (see [l4| for a list). By definition, C 
is a set of operations under which E cannot increase; 
and the bounds on Tc tell us how close the state p 
can be brought to the singlet state using only opera- 
tions in C. If C = LOCC, each p af3 can further be 

twirled, leading to the map p — ► p *+ + (1 -p)\ with 
p = {iJ-LOCc(Sobs) + 1) /3. For such states the en- 
tanglement measures can generally be computed. For 
instance, using pH |. the entanglement of formation is 
bounded by 



E f >h[- 



1 



8(1- V2)+45 0bs -5 2 , 



obs 



,2 4(72-1) 
where h is the binary entropy function 



V. CONCLUSION 

A theory of black box source characterization is a step 
towards the development of device-independent quantum 
information processing. In the present work we used only 
the CHSH inequality: already in this simple case, we 
have uncovered a rich structure, raised many problems 
and solved a few. 

In particular, the task of deriving black-box bounds 
for use in the black-box scenario in full generality is still 
open; we have been able to derive tight bounds for the 
Mayers- Yao fidelity either by restricting the dimensions 
to two qubits (JTUJ), or by restricting the state to be pure 
(fl"6]l . For arbitrary states, we do not even have a lower 
bound for the Mayers- Yao fidelity or trace distance. How- 
ever we have been able to derive unrestricted black-box 
bounds for use in other scenarios (| 19I20[) where one wants 
to ascertain how close to an ideal state it would be pos- 
sible to bring the system by local operations, possibly 
complemented by classical communication. We have also 
been able to derive unrestricted black-box lower bounds 
for the entanglement of the state. 

Our results indicate that black-box bounds put very 
stringent demands on the quality of an untrusted source, 
which could in particular have important consequences 
for self testing of quantum computers. 
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APPENDIX A: CALCULATING F MY FOR PURE 
STATES 

We begin with a state in Schmidt form 



(Al) 



while the closest state of the form |?) ® \<j>+) has Schmidt 
decomposition 



(A2) 



with (I21 = ^2i+i- For concreteness, we may assume that 
the XjS and fijS are both in decreasing order 

We first show that we may take \cj) = \a,j) and \dj) = 
\bj). Note that 

\(m\<J2 X ^k\(a J \c k )\\(b ] \d k )\ (A3) 



Let us define the matrix M by 



M jk = \(a j \c k )\\(b j \d k )\. 



(A4) 



The values | (<Xj\bk) | for various k and fixed j form a vector 
of norm 1 since \b k ) is a basis and \a,j) has norm 1. The 
same is true for the values |(6j|dfc)| and if we fix k and 
vary j instead. Thus columns (and rows) of M are formed 
by entrywise products of norm 1 vectors and the sum of 
each row and column of M is at most 1. This means that 
we can find a new matrix with positive entries such 
that M + N is doubly stochastic. Note that 



<J2 x ^k(M+ N) jk . 



(A5) 



By the Birkhoff-von Neumann theorem we may write 
M+iVasa convex combination of permutation matrices, 
thus 



M 



N = ^p m P„ 



(A6) 



with ^2 m p m = 1 and P m permutation matrices. Since 
the combination is convex, there exists some m for which 



jk 



(A7) 



The permutations merely reorder the HjS and it is easy 
to prove that the maximum is achieved when the AjS and 
fijS are both in decreasing order. Hence P m = I satisfies 
the above equation. We may achieve this by choosing the 
bases |c,-) = leu) and IcL) = 



any other bases 



so we need not consider 
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We now optimize over pj subject to the condition pii 
M2Z+1- By the Cauchy-Schwarz inequality we have 



(A8) 



with equality when p and A are collinear. Thus we set 

A22 + ^21+1 



P21 = M2Z+1 



N 



with N a normalization constant equal to 
N 



/2^(A 2 i + A 2 i + i) 



(A9) 



(A10) 



With these values, we obtain 



with LO the set of local operations that take the space 
AB to a pair of qubits. We may restrict this set to op- 
erations which only apply local unitaries and trace out 
everything but a pair of qubits to obtain 



Flo(p) > maxF{ti x {U ®Vptf ®V^),\(f> + ){<j> + \) (B2) 



where trx means tracing out everything but a pair of 
qubits. Since Fidelity only increases when a system is 
traced out we have 



F L o{p) > maxF((7 ® Vprf <8> V\ \4>){<t>\ ® \<t>+)(<t>+\) 

(B3) 

for all \<j>), and in particular for the \<f>) which maximizes 
the expression and gives Fmy(p)- Thus 



f M k(|V» = KV#}| 2 = £ 



(X 2 l + \2l + lY 



(All) 



Flo(p) > F MY (p) 



(B4) 



APPENDIX B: PROOF OF F MY < F LO WITH 
EQUALITY FOR PURE STATES 



Let p be given. Then 



F LO (p)= maxF($(p),|^ + }(0 + |) 



(Bl) 



Now suppose that p = \ip)(tp\ A B- We may write an 
operation in LO as adding a pair of ancillas and a pair 
of target qubits, applying a pair of unitaries, and tracing 
out everything but the target qubits. Thus 



F LO {Wi) = va^F{tv ABXa x b {U^V\i;)AB\W)x a x^)Y a Y b ),\<t>+){4>+\Y a Y b ) 



Applying Uhlmann's theorem, we obtain 



F £0 (|^» = max \(tp\ AB (00\x a X b (W\Y a Y b tf ®V^\cl>) ®\<t> + ) YaYb 
u,v,\4>) 



(B5) 



(B6) 



The right hand side is equal to F MY {\ip) ® |00) ® |00)) by 
definition. This in turn is equal to FmyQiP)) since the 
value of Fmy for a pure state is only dependent on the 
Schmidt decomposition, which the product state ancillas 
do no change. Thus 



F L o{\i>))=F MY {\i>)) 



(B7) 



For mixed states there exist cases with a strict inequal- 
ity. For example Fmy(j) = \, but Flo(j) = \ since the 
class LO allows us to replace the state with 1 00) . 



APPENDIX C: PROOF OF THE LOWER BOUND 

FOR F lo . 



Here we prove the bound on Flo cq. (fT9|) . The defi- 
nition of Flo i s: 

F LO = max V Tr \M k <g> N L pMl <g> N} *+ 
M k ,Nij-f L fc ' J 

M k :H A ^<C 2 ; J2 M l M k = ^A 

k 

Ni:H B ^C 2 ; ^ JV/iVj = 1 B (CI) 
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where {Mk} ({Ni}) are CP maps from Alice (Bob's) sys- 
tem to 2 dimensional spaces. In general {M^} and {Ni} 
will depend on p. 

We will explicitly describe CP map's that achieve eq. 
(fT9]) . thereby showing that it is a lower bound for Flo- 
This bound is certainly not tight. This can be seen by 
the construction we use, since the CP maps depend on 
the measurement operators A. A', B, B' , but not on 
the state p itself. Thus the CP maps do not use all the 
available information, and cannot not be optimal. 

The CP maps are constructed as follows: 

1. The operators A, A' (and B, B') are block diagonal, 
where each block is a 2 x 2 matrix. We use the 
projectors n Q <E> Hp to project onto these blocks, 
obtaining states p a p with probability p a p- The Bell 
operator in block (a, (3) has expectation S(p a p). 

2. If S(p a p) < 2, then F > 0. In this case do nothing. 

3. If S(p a j3) > 2, then carry out local rotations, such 
that after the rotations, the measurements look like 

A = cos aZ + sin aX 
A 1 = sin aZ + cos aX 



Explicitly we have: 



7T 

< - 

1 1 - 4 



(C2) 



and 



B = cos(^ + b)X + s\n(^ + b)Z 

= -^(cos b — sin b)X H — = (cos b + sin b)Z 
V 2 v 2 

B' = - cos(^ + b)X + sin(^ + b)Z 

= p(cos6 — sa\b)X H — -= (cos b + sin b)Z 

V 2 v 2 



lis; 



(C3) 



Ai(o, b) = 2vl + cos 2a cos 26 



2\f2\/ cos 2 a cos 2 6 + sin 2 a sin 2 & 



A 2 (a,6) = 2 VI — cos 2a cos 26 



(C4) 



The fidelities of the eigenvectors with the \& + ) state are 
Si = |<$ + |*i) | 2 . One finds 



/ 2 = /a = 
A + /4 = 1 

2^2 cos a cos 6 



A/ = A - / 4 = 



Ai(o, 6) 



(C5) 



The pure state on which the measurements are carried 
out can be written in the basis as: 

\f) = Cl|$l) +c 2 |$ 2 ) +c 3 |$ 3 ) +c 4 |$ 4 ) • 

The expectation of CHSH is 

(B) = S = ( I ci 1 2 - |c 4 | 2 )Ai + (|c 2 | 2 — |c 3 | 2 )A 2 

and the Fidelity is 

F=\c 1 ^h+c i ^h\ 2 

(where we take the \/7i an d y/JI to be the postive square 
roots of fi and / 4 ). 

Note that these expressions would be unchanged if we 
had a mixture of | < I > ) 2 , \&}3, and c 4 1<&) 4 +c 4 |$) 4 . Hence- 
forth we consider such a mixture. If the state is of the 
form I tp) = |$) 3 or of the form = |$) 2 , then S < 2, 
and therefore, trivially, _F > 0. 

We now concentrate on the non trivial case \ip) = 
ci|$)i + c 4 |$) 4 and 2 < S < 2^2. We will show that 
in this case 



S 
4\/2 



< F < 1 



(C6) 



The idea of the final rotations is that the operators A, A' 
(and B, B') define local bases, and we rotate the state so 
that these bases are aligned with the local bases defined 
by the state <&+. 

Note that as both the fidelity and the CHSH viola- 
tion are linear functions of the density matrix p, we can 
restrict ourselves to pure states. Furthermore, because 
of the linearity of S and F, we can focus on one block 
(a,/3). Taking the concave hull will yield the set of ac- 
cessible points. From now on we drop the indices a, (3. 

Using eqs. (|C2[ IC3|) the Bell operator takes the form 

B = \/2[cosacosb{ZZ + XX) + cosasmb(ZZ -XX) 
+ sin a cos b(ZX + XZ) + sin a sin b(-ZX + XZ)\ 

The eigenvectors of B are denoted and its eigenvec- 
tors are Ai, A 2 , — A 2 , — Ai where < A 2 < 2 < Ai < 2y2- 



Note that eq. (|C6[) and the preceding arguments give 
us the extremal points in the (S, F) plane we were search- 
ing for. Taking the concave hull yields eq. (fT9| . The con- 
cave hull can be attained by taking the angles a = b = 0, 
and as state a mixture of \cj) + ) (which has F = 1 and 
S = 2V2 and of |$ 2 ) (which has F = and S = 2). 

Proof of eq. flUgJ) . 

To prove eq. (|C6|) this recall that F = |civ / 7T+c 4 i 
and S = (|ci| 2 — |c 4 1 2 ) A 4 . We can view F = \v ■ w\ 2 
as the scalar product of two vectors v = (ci,c 4 ) and 
w = {\fj\,\fjl)- Our argument will be to fix S and to 
minimize F. 

For fixed 5, a, 6, the minimum of F is obtained when 



TA 2 



Cl = 



l + ' C4 = -Vi"^T- (C7) 



From now on we take c\ , c 4 to have this form. 
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We can then write 



F = 



2 2Ai V 2 



1 S 1 Af I S I Af 



2 2Ai V 2 



(C8) 

From now on, our aim is to choose the measurement an- 
gles a, b that minimize eq. (|C8j) for fixed 5*. 

First let us keep S and Ai fixed. Then F is minimum 
when Af is minimized. We show that this occurs when 
\a\ = \b\. 

Proof. We consider the (a, b) plane. The vector 
ft = (— cosasina(cos 2 b— sin 2 b), — cos&sin6(cos 2 a— sin 2 a)) 
is normal to the surfaces Ai = constant] and the vector 
t = (cos6sint>(cos 2 a — sin 2 a), — cosasina(cos 2 b — sin 2 b) 

is tangent to the surfaces Ai = constant. 

Recall equation (|C5[) . It then follows that t ■ 
(— sin a cos b 1 — cos a sin b) = sinasin6(sin 2 a cos 2 b — 
cos 2 a sin 2 b) is proportional to the change of Af along 
the surfaces Ai = constant. Analyzing this function, one 
finds that the minimum of Af occurs when \a\ = \b\. 
End of proof. 

We can thus replace \a\ = \b\ in eq. (|C8|) . Then when 
S > 2 the minimum of F occurs when a = b = 0. Replac- 
ing a = b = in eq. (|C8[) . this is equivalent to proving 
that when \a\ = \b\, F>± + 



Proof. We use cq. (IC8|) to rewrite the inequality 

F> I-L-fiL 



2 1 4/2 



as 



2 + ^A7" 2 VI _ iA|Vl"~-2 + i7! (C9) 

which we reorganise as 



SAf _ _S^_ 1 S 2 1 Af 2 

"2A7~4V2~ V 4 ^ 4Af V 4 4~ 



(CIO) 



Both the left hand side and the right hand side are pos- 
itive (since it is easily checked that Af/S\ > 1/2V2)- 
Hence this inequality is equivalent to its square, which 
gives: 



S 2 Af S 2 



> 1 - A/ 2 - 



v^Ai 8 
Reorganising terms yields 

2(cosa 2 -l) 2 (ff 2 -4) 
A 2 



A? 



> 



(Cll) 



(C12) 



which is manifestly true when S > 2. End of proof. 
End of proof of eq. (|C6j) . 
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